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ABSTRACT 

Magnetic reconnection at the interface between coronal holes and loops, so-called interchange recon- 
nection, can release the hotter, denser plasma from magnetically confined regions into the heliosphere, 
contributing to the formation of the highly variable slow solar wind. The interchange process is often 
thought to develop at the apex of streamers or pseudo-streamers, near Y and A-type neutral points, 
but slow streams with loop composition have been recently observed along fanlike open field lines ad- 
jacent to closed regions, far from the apex. However, coronal heating models, with magnetic field lines 
shuffled by convective motions, show that reconnection can occur continuously in unipolar magnetic 
field regions with no neutral points: photospheric motions induce a magnetohydrodynamic turbulent 
cascade in the coronal field that creates the necessary small scales, where a sheared magnetic field 
component orthogonal to the strong axial field is created locally and can reconnect. We propose that 
a similar mechanism operates near and around boundaries between open and closed regions inducing 
a continual stochastic rearrangement of connectivity. We examine a reduced magnetohydrodynamic 
model of a simplified interface region between open and closed corona threaded by a strong unipolar 
magnetic field. This boundary is not stationary, becomes fractal, and field lines change connectivity 
continuously, becoming alternatively open and closed. This model suggests that slow wind may origi- 
nate everywhere along loop-coronal hole boundary regions, and can account naturally and simply for 
outflows at and adjacent to such boundaries and for the observed diffusion of slow wind around the 
hcliospheric current sheet. 

Subject headings: Magnetic reconnection — magnetohydrodynamics — solar wind — Sun: corona — 
Sun: magnetic topology — Turbulence 



1. INTRODUCTION 

A topic of recent interest is magnetic reconnection 
between open and closed field lines at the interface 
between coronal holes and loops, dubbed "interchange 
reconnection" (Figure [TJ IR hereafter). This mech- 
anism can contribute mass, heat, and momentum to 
the solar wind, with numerous heliqspheric implications 
(iFisk et al.lll999HFisk fc SchwadrorfeOOlHCrooker et al.l 
2002t iDahlburg fc Einaudil 120031: lAntiochos et all 120071: 
Owens et al.l 12008: Ed mondson et al.l l2009t iLinker et al.l 



2011; Titov et al. 2011; Masson et al. 2012). 



The solar wind may be classified as "fast", when ve- 
locities v exceed, say, 600 kms -1 , and "slow", when 
v < 500 kms -1 . The steadier fast wind originates in po- 
lar coronal holes (dark X-ray regions) and similar open- 
field regions closer to the equator, pro pagating radially 
into interplanetary space (|Zirkerlfl977l ). The more spa- 
tially and temporally intermittent slow wi nd originate s 
in and around the coronal streamer belt (jWangl Il994f l . 
where a significant population of closed-field structures 
is found. The heliospheric current sheet (HCS) is al- 
ways embedded within slow wind, which surrounds it in 
a region spannin g about 30° in latitude near solar mini- 
mum conditions (iGosling et al.lfT98Tt iBorrini "ellllllQSlI : 
iWinterhalter et al.lll9941 ). 

Fast and slow wind differ in their plasma composi- 
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tion. Generally the fast wind composition is similar 
to that of the photosphere, while slow wind compo- 
sition is similar to that of coronal loops, with com- 
parable abundances ratios of low to high first ion- 
ization potential (FIP) elements and ions with differ- 
ent charge states (e.g Q 7+ /0 6+ ) (IGeiss et all 119951: 
iZurbuchen eTaTl [2002t iFeldman fc WidinsJ 120031) . IR 
allows field-line connectivity to change from closed to 
open, thus releasing coronal loops plasma into the he- 
liosphere. This has been suggested as a pri mary mech- 
anism for the formation of the slow wind (jWang et al.l 
Il99l lAntiochos et aLllMtl . 

Recent observations motivate to advance our under- 
standing of the physical processes at the root of IR. 
Outflows along open fanlike field lines, at the edges 
of active (closed) regions, have been observed by Hin- 
ode EUV imagin g spectrometer (EIS) me asurements 
(jSakao et al.ll2007| ). iBrooks fc Warrenl (|2011| ) found that 
the composition of these outflows is that of coronal loops 
and established a link with slow wind detected in situ 
by the Solar Wind Ion Composition Spectrometer on 
board the Advanced Composition Explorer (ACE). Out- 
flows in the streamer belt region ar e observed by the 
LASCO-C2 white light coronagraoh ilWang et al.l 12011 
and STEREO imagers ([Howard et al.ll2012D . suggesting 
that mixing and dynamics contribute to the average ob- 
served configuration. 

In most prior work, IR has been thought to occur only 
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Figure 1. Cartoon of interchange reconnection in a turbulent 
Corona. Reconnection takes place in point X, along the bound- 
ary between open and closed corona, with closed field line A-X-B 
opening into A-X-C, while the field line traced from point B closes 
at a different location, point D, forming the closed field line D- 
X-B. Swinging field lines depict magnetic fluctuations (pictorially 
exaggerated). 

in special topological locations, at the apex of streamers 
and pseudo-stream ers corresponding to Y or X-points 
ijWang et al.ll2012l ). where field lines of opposite polarity 
c an reconnect in a n eutral point with B = 0. 

I Wang et al.l (|1998f ) proposed that convective field line 
shuffling can trigger IR around the cusp region, resulting 
in the outward propagation of density enhanced blobs. 
But the physical mechanism leading to and allowing IR 
remains undetermined. In fact, the small value of re- 
sistivity in the solar corona implies that magnetic field 
lines are frozen in the plasma except where very small 
scales are present, i.e., strong currents with an apt lo- 
cal magnetic field topology for magnetic reconnection to 
occur. 

Numerical simulations (lEinaudi et al.l 119961 : 



Dmitruk & Gomez 1997; Rapp azzo et al 



Parker model for coronal heating (IParkerl ll972Tfl988h 



20071) o f 



have shown that the continuous shuffling of magnetic 
field lines' footpoints by photospheric convective mo- 
tions induces a magnetohydrodynamic (MHD) turbulent 
cascade in the unipolar closed coronal field with no null 
points. This cascade transfers energy from the large 
to the small scales, driving field-aligned current sheets 
that are continuously formed and dissipated, where the 
magnetic field component orthogonal to the strong axial 
field is sheared, i.e., its field lines are locally oppositely 
directed, and can reconnect (nanofiares). 

Therefore the unipolar closed field lines of coronal 
loops continuously change connectivity due to this dy- 
namical activity. Furthermore, in view of ubiquitous 
presence of magnetic fluctuations in this scenario, at 
each instant of time the magnetic field lines admit a 
random character due to field line random-walk (F LRW: 
Uokipii fc Pa7kei1[l968l fl969t iMatthaeus et al.|[2007T) . In 
this environment, magnetic connectivity is very complex 
and changing. 

Here we suggest that similar dynamics take place ev- 
erywhere at the boundary between open and closed re- 
gions where turbulent IR can occur stochastically (Fig- 
ures HE]), naturally accounting for the observed flows 
along and around these boundaries, includ ing those at 
adjacent active region edges observed by iSakao et al.l 




Figure 2. Snapshot of the simulated interface between closed 
(straightened out loop) and open regions of the solar Corona. 
Field lines are line-tied at the photospheric plane 2 = 0, where 
convection-mimicking motions (shown in contours) are applied at 
x < 1 and vanish for x > 1. Closed field lines return to the plate 
z = 10 for x < 1 where they are line-tied to a motionless photo- 
sphere, while for x > 1 an open boundary is realized. The plane 
x = 1 is the original boundary magnetic surface between open and 
closed regions at t = 0. 

(|2007f ). that cannot be explained by IR at the streamer 
apex. 

In this paper we investigate the dynamics of IR at the 
interface between open and closed corona, with photo- 
spheric convective motions shuffling the magnetic field 
lines' footpoints. For a simple first demonstration, we ap- 
ply photospheric motions only to the (originally) closed 
region, so that no waves or turbulent dynamics are ex- 
cited directly by photospheric motions along the origi- 
nally open field lines. 

2. MODEL AND GOVERNING EQUATIONS 

We model the interface region in Cartesian geometry 
(Figure [2]), with a straightened loop juxtaposed with an 
open-field region. Curvature effects are neglected. The 
computational box spans 0<a;<2, 0<y<l, and 
< z < 10 (lengths are normalized by £* = 60 Mm, e.g., 
z spans 600 Mm). The plane z = represents the pho- 
tosphere where both open and closed field lines are line- 
tied. On the z = plane, in the x < 1 region, convective 
super-granular motions are modeled by imposing a large- 
scale velocity pattern with all modes with wave-numbers 
between 3 and 4 excited with random amplitudes and 
normalized to have an r.m.s. value ~ 0.5 kms -1 , that in 
physical space correspond to distorted vortical stream- 
lines with length £ c ~ 15 M m. This model of foot point 
motion is similar to that of Rappazzo et~aT1 (pfM ). and 
is illustrated in the contours at the bottom of Figure [2j 
On the remaining region of the z = plane, where x > 1, 
the velocity vanishes. The upper plane at z = 10, in the 
region x < 1, represents the photospheric plate where 
closed loop field lines return to, and are line-tied to a 
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motionless photosphere. On the section of the z = 10 
plane having x > 1 , an open boundary i s realized, impos- 
ing non-reflecting boundary c onditions (|Thompsonlll987l 
[19901 rVanajakshi et all [1981 . i.e., wave-like signals are 
allowed to propagate out toward z > 10 with n o reflec- 
tion toward z < 10 (e.g.. iRappazzo et al.l 12005b . Along 
x and y periodic boundary conditions are used. 

The above system is threaded by a strong and uniform 
unipolar magnetic field Bo = Bo e z along z. The field 
lines traced from the bottom photospheric plane z = 
are considered to be either closed when they map to 
the top z = 10 plate with x < 1, or open for x > 1. 
Because of a large assumed conductivity and line- tying, a 
field line must undergo magnetic reconnection to change 
connectivity. Field lines traced from the plate z = 10 
with x < 1 map the actual closed region in the z = 
plane. Likewise those traced from z = 10 with x > 1 
map the open region back to the z — plane. 

As in previous work the dynamics are integrated 
with the (nondimensional) equations of reduced magne- 
tohydrodynamics fMHD) ([ K adomtsev fc Pogutsd 119741: 
Strauss 1976; Montgomery 1982), well suited for a plasma 
embedded in a strong axial magnetic field: 

ftu + u Vu^ -VP + b ■ Vb + c A d z b + ^— '- V 2n u, (1) 

ne n 

<9tb + u- Vb = b Vu + c,9 z u+^— '- V 2n b, (2) 

Re n 

with V • u = V • b = 0. Here, gradient and Lapla- 
cian operators have only transverse (x-y) components as 
do velocity and magnetic field vectors (u z = b z = 0), 
while P is the total (plasma plus magnetic) pressure, c A 
is the Alfven velocity of the axial field (B / 'iirpo), the 
plasma is assumed to have uniform density po- In the 
simulation presented here c A = 200 (velocities are nor- 
malized by u* = 0.5 kms" 1 ), and the numerical grid has 
512 x 256 x 120 points to achieve the long duration of 
~ 1,800 Alfven crossing times ta = L z /c A (L z = 10 
is the loop length), corresponding to ~ 180 nonlinear 
times, necessary to acquire significant statistics. We use 
hyperdiffusion with n = 4 and R4 ~ 5 x 10 16 that elim- 
inates diffusion at the large scales, a critical feature in 
this kind of simulations that otherwise reach a diffusive 
regime (see, als o for a more detailed description of the 
numerical code, IRappazzo et al.ll2008[ ). 

3. RESULTS 

Initially photospheric motions induce magnetic fluctu- 
ations in the closed regions that grow linearly in time. 
A nonlinear, turbulent sta ge is attained around time 
t nl ~ 10 ta (for details see IRappazzo et al. 2008). The 
dynamics subsequently leads to many field-aligned cur- 
rent sheets where magnetic reconnection occurs. These 
structures are highly dynamic, crossing a transverse cor- 
relation length (here, the super-granulation scale £ c ) in 
approximately a nonlinear time-scale t nl . In this dy- 
namic sea of current sheets, roughly half of those forming 
within a correlation length from the open-closed bound- 
ary will encounter it, thus changing its location and in- 
ducing changes in field-line connectivity. 

Figure [5] shows field lines at t ~ 163.59 ta originat- 
ing at z — near the initial open-closed boundary at 
t = (x = 1), but which cross the boundary prior to 
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Figure 3. Field lines' footpoints in the plane z = 10. Field lines 
are traced from points with x = xo in the plane z = 0. Their 
footpoints' orthogonal displacement relative to their origination 
point (central circle) is shown for four selected values of xo. The 
boundary between open and closed regions is at x = 1. 

arrival at z — 10. This opening and closing of field lines 
is caused by the highly dynamic current sheets and re- 
connection described above. This is bursty and stochas- 
tic, and increasingly so with higher Reynolds numbers 
(jServidio et al.l[2W9l ). 

To quantitatively understand the impact of this kind 
of IR we analyze the statistical properties of field lines. 
We ask what is the fraction of time spent in open/closed 
regions, or the probability for a field line traced from 
a point x = xq from the boundary in the photospheric 
plane z = to be closed or open when it arrives at the 
top plane z = 10. The latter condition corresponds to 
X < 1 or x > 1 respectively. 

To address this, we trace field lines from points in the 
plane z = with x = xq fixed. They are traced at 320 
different times, separated by At = 5 ta, corresponding to 
approximately half a nonlinear time. To increase statis- 
tics the field lines are computed in 40 equally spaced 
points along y (along this direction points are statisti- 
cally equivalent). Figure [3] shows the footpoints' dis- 
placement in the x-y plane for 12,800 field lines in the 
plane z = 10. Four cases correspond to four selected 
initial distanc es xp from the bo undary. The field line 
tracing code (jDalena et al.l [20121) employs a fifth order 
Runge-Kutta with adaptive step-size, and second order 
interpolation. 

Field lines traced from the middle of the originally 
closed and open regions, i.e., Xq = 0.5 and 1.5, exhibit an 
isotropic distribution of footpoints of different extension 
because in the closed region magnetic fluctuations are 
stronger. No waves are injected from the region z — 
with x > 1, but magnetic field fluctuations "leak" from 
the closed region, where magnetic forces push the mag- 
netic islands against each other and these forces are un- 



4 



Rappazzo et al. 




y t/t A = 168.69 



0.6 



0.8 



1.0 



1.2 



1.4 



Figure 4. Probability density functions of the footpoints shown 
in Figure[3]for xo = 0.9 and xo = 1.1, across the boundary between 
open and closed regions at x = 1. 

balanced by weaker fields in the open region, with the 
energy density reducing across the original boundary at 
x = 1. This drop in turbulence intensity along x is the 
cause of anisotropy for the footpoints distributions for 
Xo = 0.9 and 1.1 that are narrower along x. 

Field lines traced from xq — 0.5 and 1.5, more distant 
from the initial open-closed boundary, never change con- 
nectivity from open to closed (or vice- versa). But those 
traced from xq = 0.9 and 1.1 do, and their p.d.f.s are 
shown in Figure |U They have a probability of ~ 20% to 
change connectivity from their initial type at time t = 0. 
From the dynamical point of view the probability can 
be seen as the fraction of time the field line is closed or 
open. 

Given the complex character of the magnetic field b, 
which becomes broad band in space, while also evolving 
in time, it is natural to try to describe the spreading (and 
changes of connectivity) as a diffusive process. In fact 
there are two related diffusive processes at work. First, 
in Figure [2] we see that at a fixed instant of time, field 
lines wander randomly due to fluctuations. Second, for 
the spreading found in the above numerical experiment 
shown in Figures [3] and IU the time dependence of the 
turbulence contributes an additional randomizing effect. 
The diffusive nature of the field line spreading is verified 
by the computation (not shown) of the mean square dis- 
placement (Ax 2 ), which increases linearly with height in 
the manner of diffusion: 



(Ax 2 



2Dz. 



(3) 



We found this linear scaling for all values of xo, with 
diffusion coefficients, e.g., D = 93.6, 56.82 and 14.46km 
for Xo = 0.5, 1, and 1.5 respectively (normalized to £* = 
60 Mm). 

This type of randomization of magnetic field lines is not 
often associated with coronal fields (usually modeled as 
smooth and laminar), where field lines may be line-tied at 
both ends. This diffusive rearrangement of connectivity 
requires magnetic reconnection to occur. A quantitative 




Figure 5. Open-closed regions boundary magnetic surface. 

theory for this space-time diffusion of field lines appears 
to be tractable but lengthy, and we will address it in a 
subsequent paper (Ruffolo et al., in preparation). 

Two features of the diffusion theory are pertinent at 
present. First, the expected (Ax 2 ) at height z due to 
single time randomization - the FLRW effect - and the 
expected additional mean square displacement at height 
z due to the time dependent changing of field lines, are 
of the same order. Second, there are two broad classes 
of FLRW diffusion coeffi cients, say D a i ~ X z (b /B ) 2 , 
the quasi-linear re sult dJokipii k, Parker! 11968ft . and 
D B ohm ~ X ± b/B a (|Ghilea et all 1201 If ). Here X z ~ 5 
and A_l ~ 1/4 are suitable coherence scales in directions 
parallel and perpendicular to Bo, with bj Bq ~ 2% in 
the closed region. Note that the present case of reduced 
MHD requires the ordering bX z /(B X±) ~ tm/t± ~ 1, 
from which DBohm ~ D q i. Consequently we anticipate 
that the observed diffusive spreading (Figures [2][4j) is 
characterized by a diffusion coefficient on the order of 
the quasi-linear result. 

Also important is the evolving structure of the bound- 
ary between open and closed regions. At time t = the 
boundary is simply the plane x = 1 (Figure [2]), but at 
later times the magnetic surface has to be computed. 

In reduced MHD, with uniform axial field Bo = £>oe 2 
and transverse fluctuations b, the magnetic surface co- 
ordinate ifi obeys the magnetic differential equation 



d z i> = -— b • V?/>, 
Bq 



(4) 



where the right side involves only the components of 
V transverse to Bo, with initial condition tp(x,y,z = 
10) = sin(-Kx). The numerical code employs a third or- 
der Runge-Kutta, quadratic interpolation, and adaptive 
step-size. 

Like a passive scalar, sol utions to equation (|4"|) c an ac- 
quire a complex structure ((Matthac us et al] [l995V The 
boundary magnetic surface ip — separates the two topo- 
logically different regions, and in Figure [5] it is shown at 
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time t — 168.69^. Its structure is fractal, appears like 
a pleated drape with many intricate folds, but for a con- 
tinuous held b it does not tear no matter how folded 
it is, although numerically a small diffusivity removes 
the smallest-scale folds and minor tearing occurs. The 
boundary magnetic surface evolves in time and on the 
average its map on the plane z = has an excursion in 
x given by twice (Ax 2 ) 1 / 2 = {2DL Z ) 1 / 2 (equation ©), 
where L z is the loop length. 

4. CONCLUSIONS AND DISCUSSION 

Previous studies have shown that IR is required, e.g., 
to explain the quasi-rigid rotation of coronal holes in 
presenc e of the underlying photospheric differential ro - 
tation (IWang et all [1991 ILionello et all 120051 [20061) 
but this approach (se e also iSchriiver fc De Rosal f2003: 
IWang fc Sheelevl l200l assumes a quasi-steady coronal 
response to photospheric evolution. In the prevailing 
view that coronal interchange occurs at the apex of 
streamers an d pseudo-streamer s in correspondence of Y 
or X-points (Wang et al. 2012), all previous simulations 
and modeling have used smooth large-scale fields that 
c ontain neutra l point s (B = 0). 

IWang et al.l (|1998l ) anticipated that held line foot- 
points shuffling may promote IR at the boundary be- 
tween open and closed regions. The present model pro- 
vides a specific mechanism for this to occur, modeling 
IR as component reconnection that may occur all along 
the magnetic interface between coronal hole and loop 
threaded by a unipolar held that contain no true neutral 
points, a nd extends t h e ran ge of occurrence of IR pro- 
posed by IWang et al.l (|1998l ). These reconnection sites 
are well known in the context of nanohare models, but 
their role in interchange has not been previously empha- 
sized. Here we simulate only a small volume of the open- 
closed region interface to employ a higher spatial resolu- 
tion. This allows the development of MHD turbulence, 
with its associated magnetic fluctuations in the coronal 
held (Figures[T]and[2]), naturally induced by photospheric 
motions shuffling the held lines' footpoints. 

Turbulent IR renders the boundary between open and 
closed regions dynamic (Figures [3][5|). The boundary 
fluctuates continuously with an average displacement of 
the order of the super-granulation scale £ c ~ 15 Mm. 
IR can then inject loop plasma along the boundary 
and in the fanlike regions adjacent to closed regions, 
where slow strea ms with loop composition have been 
recen tly observed (jSakao et al.l [2007: Br ooks fc Warrenl 
1201 If) , providing an alternate mechanism to account for 
the plasma composit ion at the edges of active regions 
(|Cranmer et al.ll2007f) . and additional momentum , mass 
and e nergy for the streams originating from there (jWand 
fl99l . 

In a realistic geometry the field lines originating from 
this small fanlike regions expand super-radially in the 
heliosphere and map in an extended region around the 
HCS. Thus flows due to turbulent IR naturally diffuse 
away from the HCS, overcomin g the restriction of the 
smooth field model proposed bv lAntiochos et afl (|2011f l 
that admit diffusion only for streams originating from 
narrow open flux channels connecting two coronal holes. 

In summary, when held lines' footpoints are shuffled 
by photospheric motions, component magnetic reconnec- 
tion is expected to occur in unipolar loop and open field 



regions and near the boundaries between them. This 
stochastic IR is likely to operate all along these bound- 
aries and adjacent regions, where closed and open field 
lines can then continuously change connectivity. On this 
basis we suggest that plasma and energy transport along 
these magnetic held lines may be an important factor in 
generating the slow wind, and in broadening the regions 
in which compositional and other properties are mixed 
in the solar wind. 

In the future we plan to extend this work to more re- 
alistic reduced and full MHD models that include curva- 
ture and expansion effects and alternate boundary con- 
ditions, allowing us to determine the relative importance 
of apex neutral point IR and stochastic component IR. 
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